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ABSTRACT 

Whilst the Vela pulsar and its associated nebula are often considered as the 
archetype of a system powered by a ~10^ year old isolated neutron star, many 
features of the spectral energy distribution of this pulsar wind nebula are both 
puzzling and unusual. Here we develop a model that for the first time relates 
the main structures in the system, the extended radio nebula (ERN) and the X- 
ray cocoon through continuous injection of particles with a fixed spectral shape. 
We argue that diffusive escape of particles from the ERN can explain the steep 
Fermi-LAT spectrum. In this scenario Vela X should produce a distinct fea- 
ture in the locally-measured cosmic ray electron spectrum at very high energies. 
This prediction can be tested in the future using the Cherenkov Telescope Array 
(CTA). If particles are indeed released early in the evolution of PWNe and can 
avoid severe adiabatic losses, PWN provide a natural explanation for the rising 
positron fraction in the local CR spectrum. 

Subject headings: Gamma rays: general. Gamma rays: ISM, (ISM:) cosmic rays, 
(Stars:) pulsars: individual: PSRB0833— 45 



1. Introduction 



The proximity (distance: 290 pc, measured by parallax Caraveo et al. 2001 Dodson 



et al. 2003) of the Vela pulsar (PSRB0833— 45) has allowed it to be studied in great detail 
across the whole electromagnetic spectrum. PSRB0833— 45 (spin-down power E = 7 x 
]^g36 erg/s, characteristic age Tc = 11,000 years (Manchester et al. 2005)) is embedded in 



a diffuse radio nebula ( 


Large et al. 


1968 


). Higher an| 


1991 


Duncan et al.||1996 


) and X-ray 


(Kahn et al.||1985; 



gular-resolution radio 1 


Dwarakanath 


Aschenbach et al.||1995 


) observations 
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of the region established the presence of (1) a large circular shell (about 8° in diameter) - 
the Vela SNR - and (2) an extended radio nebula (ERN) of size 2° x 3° in the centre - 
dubbed Vela X (see Figure [l] left). Based on the high degree of polarisation and the fiat 



radio spectral index Weiler & Panagia (1980) suggested that Vela X is a pulsar wind nebula 



(PWN) powered by PSRB0833— 45. Such PWNe act as calorimeter, containing the time- 
integrated particle outflow and magnetic field from the pulsar. The radio emission is far 



from uniform, exhibiting a network of filaments throughout the nebula (Frail et al. 1997) . 



The brightest of these non-thermal filamentary structures in Vela X emerges to the south 



of PSRB0833-45 and has a size of 45' x 12' (Figure [T 



this bright radio filament was detected using ROSAT (Markwardt & Ogelman 1995) (see 



right). The X-ray counterpart to 



Figure [T] right) and dubbed "the cocoon". It is thought to be the result of the interaction of 



the PWN with the SNR reverse shock (Blondin et al. 2001). 



X-ray observations performed by Mangano et al. (2005) found that within 1' of the 



Vela pulsar the 3-10 keV spectrum softens with increasing distance, a signature of electrons 
cooling via synchrotron emission. BeppoSax data demonstrated that the nebula emission 
in the inner 12' in radius extends up to ~200 keV (Mangano et al. 2005) with most of the 



emission coming from within 4' of the pulsar, suggesting particle acceleration close to the 
pulsar (black circle in Figure [T| . 

H.E.S.S. detected TeV 7-ray emission from Vela X exhibiting an atypical spectrum with 
energy fiux peaking at ~10 TeV QAharonian eraL][2006| ) . The TeV emission is dominated by 
an elliptical region of 58' x 43', extending significantly beyond the X-ray cocoon, yet smaller 
than the ERN (see Figure [l] right). The latest H.E.S.S. results (Dubois et al. 2009), show 
significant, but low-surface brightness TeV emission beyond this region. Perhaps surprisingly, 
the 7-ray spectrum in the outer parts is identical to that of the central region. The H.E.S.S. 
detection of emission from the cocoon provided the first robust estimate of the magnetic 
field strength in this region: ~ 4/iG. No excess TeV emission was seen from the hard X-ray 
emitting region close to the pulsar, implying a magnetic field of -B > 50/iG in the immediate 
vicinity of the pulsar, consistent with the cooling signature detected with BeppoSax. 

Finally, Fermi-LAT observations in the Vela pulsar off-pulsa data established the pres- 
ence of an extended 7-ray emitting structure at energies above 800 MeV spatially coincident 
with the ERN (see Figure |l| left) with a spectral index of ~ 2.4 ( |Abdo et al^|2010D . The 



two-peaked nature of the 7-ray spectral energy distribution (SED) led the authors to con- 



clude (following an earlier suggestion by de Jager et al. 2008) that there are two distinct 
populations of electrons, one responsible for the radio and GeV 7-ray emission and the other 
for the X-ray and TeV emission. The magnetic field derived through these observations is 
~ 4/iG. 
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Fig. 1. — Left: Radio image of a 3° x 3° region around Vela (black triangle), from 8 GHz VLA 



data (Frail et al. 1997). Overlaid are Fermi-LAT 4,5 and 6 a significance contours (green) 



above 800 MeV (AbdoetaL 2010). Right: ROSAT (0.9-2.4 keV) count map (Markwardt fc 



Ogelman 1997), smoothed with a Gaussian kernel of 2.25', saturated at 10% of the pulsar 



peak intensity. Overlaid are H.E.S.S. smoothed excess contours (25, 50, 75% of the peak 
emission) ( Aharonian et al.]|2006 ). The black circle indicates the 4' region where most of the 



hard X-ray emission originates (Mangano et al. 2005). 



Despite the wealth of experimental data and several theoretical studies, many questions 
remain about the relationship between the emission seen on different scales in Vela X. In 
particular, there are several important observational facts for which no satisfactory explana- 
tion has yet been put forward. Firstly, the steep GeV spectrum, indicating an unexplained 
absence of > 100 GeV energy particles from the ERN. Secondly, the relative dimness of the 
TeV Nebula. Comparing the spin-down luminosity of the Vela pulsar to the TeV luminosity 
above 500 GeV (10'^^ erg/s) yields an efficiency in the cocoon of 0.01%, significantly lower 

. In addition, the lack of spectral variations 
across Vela X above 1 TeV is not readily explained in any model in which the TeV size is 
limited by cooling. 

Here we show that a scenario in which the high energy particles have diffused out of 
the ERN, and where the cocoon is a relatively recent feature of the PWN, may answer these 
questions. This scenario has implications for the locally measured electron spectrum at high 



than typical TeV PWNe (Aharonian et al.||2006 
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energies. It has been pointed out (Shen 1970) that only very local sources can significantly 
contribute to the fiux of cosmic ray electrons above 1 TeV, with Vela (and Geminga) long 
seen as likely candidates (see e.g. Nishimura et al.||1980 Atoyan et al.||1995{ Kobayashi et al. 



2004). The recent measurements of a rise in the positron fraction above ~ 10 GeV (Adriani 



et al. 2009 Ackermann et al. 2011) indicate that a relatively local source of positrons must 



be present. Pulsars and their nebulae as plausible candidates (see e.g. Chi et al. 1996 



Biisching et al. 2008 Profumo 2008 Yiiksel et al. 2009), given that we see high-energy 



electrons/positrons radiation at TeV energies. However, the details of the escape of particles 
from the PWN region have been left open. Here we demonstrate for the first time how the 
signature of escape may be directly visible in the SED of Vela X. 



2. Model 

We explore the possibility that particle propagation effects dominate the spectral and 
morphological features of Vela X. We attempt to reconcile the properties of the TeV and GeV 
nebula by considering an electron injection from the pulsar wind termination shock which 
represents a fixed fraction of the spin-down power of the pulsar and has a fixed spectral 
shape. The two 7-ray-emitting zones of Vela X are discussed in turn. 



The TeV Nebula: Whilst the TeV emission of Vela X is now known to extend over 
much of the ERN region, the bulk of the TeV-emitting particles appear to be confined to 
a small fraction of its the volume. This emission is peaked in the cocoon, but a consistent 
spectral shape is measured throughout the ERN. This suggests that the electron population 
of the cocoon is not cooling-limited in extent and that the probability of escape from the 
system is very small and/or energy- independent. These properties appear consistent with the 
paradigm of (at least young) PWNe as internally advection-dominated systems from which 
particles cannot easily escape dAharonian etaLlpg?! ). The energetic particle content should 
then be consistent with the power supplied by the pulsar during the period of injection into 
the cocoon. As the energetic particle content of Vela X is well-determined by TeV and X-ray 
observations, the age of the TeV cocoon can be expressed as: tcoc = -Ecoc/(-^e) ~ 70yrs/e, 
with Ecoc ~ 1.5 X 10"^^ erg (Abdo et al. 2010). Whilst there are considerable uncertainties 
on the efficiency factor e of converting spin-down power into relativistic electrons for any 
plausible value of e, the cocoon is much younger than the system as a whole. This is 
consistent with an interpretation of the cocoon as a feature arising from the interaction with 
the SNR reverse shock (Blondin et al. 2001). The maximum energy of the cocoon particles 
is constrained by the H.E.S.S. measurements to ~70 TeV (Aharonian et al. 2006). This 
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maximum energy has been interpreted as a consequence of cooling of higher energy particles 
in the cocoon (de Jager 2008 Abdo et al. 2010). However, the synchrotron lifetime for 
70 TeV electrons is tgync = 10'^(-B/4/iG)~^ years, consistent with the age of the Vela pulsar 
but not with the lifetime of the cocoon. We hypothesize that the maximum energy of cocoon 
particles is determined by cooling in the high-S-field (~ lOO/iG) region in the immediate 
(~ 0.5 pc) vicinity of the pulsar before injection into the low--B-field cocoon region. The 
O(lOyr) residence time required for a cooling break to appear at 70 TeV for particles emerging 
from the high B-field region is consistent with expectations for the post-shock flow (see e.g. 



Mangano et al. 2005). 



The Extended Radio Nebula It seems likely that the ERN contains particles injected 
over the lifetime of the system, up to the evolution stage at which the cocoon appeared. In 
contrast to the cocoon, evolutionary effects in the spin-down power of the pulsar become 
critical in determining the available energy. The total energy injected into the ERN depends 
critically on the assumed birth-period of the pulsar (Po)- The energy in relativistic particles 
is ~ 2 X 10^^e(Po/30ms)~^. The electron content of the ERN has been estimated to be 
5 X 10^^ erg (Abdo et al. 2010) using the similarity in radio and GeV 7-rays. Birth-periods 
as high as Pq = 40 ms have been discussed (van der Swaluw & Wu 2001) to account for the 
ratio of PWN radius to SNR radius but require e close to 1. No explanation has been put 
forward for the absence of > 100 GeV particles in the ERN required by the steep Fermi-LAT 
spectrum. Given the current-day magnetic field of ~ 4/iG, the cooling times for particles 
emitting in the GeV-range are very large and synchrotron cooling seems implausible. While 
strong evolution of the magnetic field in PWNe is possible (see e.g. Gelfand et al. 2009), 
an implausibly high effective time-averaged B-field of > lOO/iG within the whole ERN is 
required to explain the Fermi-LAT GeV spectrum with a cooling break. Since there is no 
obvious theoretical reason why the maximum accelerated energy from the pulsar should be 
time- dependent, the remaining explanation is that the bulk of the > 100 GeV particles have 
escaped from the ERN into the interstellar medium (ISM). Whilst confinement of particles 
in PWNe is thought to be effective in their early evolution it is likely that the interaction 
with the SNR reverse shock which seems to have appeared in Vela X a few thousand years 
ago (Blondin et al. 2001; Gelfand et al. 2009) brings an end to the effective confinement 
of electrons. Such a reverse shock interaction is expected to disrupt the PWN sufficiently, 
through e.g. the growth of Rayleigh- Taylor instabilities, that diffusion of particles out of 
the PWN becomes possible. This explanation for the absence of > 100 GeV particles in the 
ERN requires a diffusion coefficient of D = i?p^/2tkyr ~ xlO^^ cm^ s~^{R in parsec, t in units 
of kyears). For typical parameters in Vela X this is of order 1000 times slower than that 
inferred for cosmic rays in the ISM, but much faster than Bohm diffusion for the magnetic 
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field in the ERN (|Aharonian||2004D 



Calculation For the calculation of the present-day SED of Vela X and the energy and time- 
dependent fiux of escaping particles we employ a Monte-Carlo approach where propagation 
and energy losses of individual particles are considered. Following earlier studies we neglect 
adiabatic losses due to the likely return of most of the lost energy at the time of reverse 



shock crushing of the PWN (see e.g. de Jager et al. 2008). Synchrotron and Klein-Nishina 
(KN) inverse Compton (IC) losses are calculated at each time step. For the calculation of IC 



cooling/radiation we adopt the soft-photon fields used by de Jager et al. (2008) and Abdo 



et al. (2010), with black-body spectra of temperatures 25 K and 6500 K and energy densities 
of 0.4 eV cm~^. In the first stage of evolution of the system confinement is assumed to be 



100% effective and a one-zone approach (following the treatment in Hinton & Aharonian 



2007) is employed. A canonical braking index of 3 (de Jager 2007) is assumed over the 



evolution of the pulsar, rather than the value measured for the recent past (determined to 



be 1.6 ± 0.1 by Dodson et al. 2003). The choice of this value has very little impact on the 
final model curves. For Pq = 30ms, the implied age of the system is tgys = 10 kyears. At 



tr 



70%tsys (the age when the reverse shock is thought to start interaction with the 



ERN, see e.g. Blondin et al. 2001), diffusive transport is switched on, with particles initially 
uniformly distributed in a spherical PWN, with continued injection until tsys— ^coc (^coc = 230 
years, e = 0.3), from which time particles injected at the termination shock are assumed to 
be confined within the cocoon. This is consistent with simulations (e.g 



van der Swaluw et al. 



2004) in that when the reverse shock has collided with the PWN, the pulsar no longer injects 



particles into the "relic PWN" but forms a new PWN. Diffusion is assumed to be a fixed 
factor S faster than Bohm diffusion, with the magnetic field in the nebula falling continuously 
with t~^-^ at late times. The free parameters of the model are the pulsar birth period Pq, 
the efficiency e of conversion of spin-down power to particles in the MeV-100 TeV domain, 
the factor S and the current-day ERN magnetic field. For the TeV nebula the magnetic field 
strength and cocoon age form the additional free parameters. To model propagation and 
energy losses of the electrons released into the ISM we adopt an analytical solution of the 



diffusion-loss equation (Atoyan et al. 1995), assuming synchrotron and IC radiation as the 



dominant loss mechanisms. Synchrotron losses assume a local ISM magnetic field strength 
of 5 fiG. The IC losses were calculated using the full KN description, assuming 4 radiation 
fields (CMB, IR, G-K stars and 0-B stars) with a total radiation density of 1 eV cm~^. 



- 7- 



3. Discussion 

Figure [2] shows measured and calculated SEDs for the ERN (red) and the TeV cocoon 



(blue). There is a striking inconsistency between the WMAP measurement shown in (Abdo 



et al.|[20l0 ) and earlier radio measurements (Alvarez et al. 2001 Hales et al. 2004). We 



attempt to approximately match the spectral shapes, rather than the ~ 10"'^ eV (50 GHz) 
normalisation, given the apparent factor ~ 3 disagreement. Data from Planck should resolve 
this issue in the near future. The main impact on our model is a change in the current-day 
magnetic field. 

As previously discussed, cooling effects are not important inside the cocoon (blue 
curves). The spectrum in the cocoon is therefore unmodified from the injection spectrum. 
The assumption is that the cooling of the particles inside the cocoon has happened close 
to the pulsar, consistent with the low magnetic field (~ 4/iG) in the cocoon and the much 
higher field (~ lOO/xG) close to the pulsar. Thus, this curve is similar to the time-independent 



one-zone IC model shown in Aharonian et al. (2006). 



The expected ERN flux (red) in the absence of particle escape for constant efficiency 
e is shown as dashed lines in Fig. |2| clearly demonstrating how dramatically different the 
observed spectrum in the Fermi-LAT energy range is. In this case the peak energy and 
fluxes are determined by particle cooling. The solid lines show the energy-escape based 
models under the approximation of spherical symmetry and a 7 pc radius. Different curves 
illustrate the effects of different diffusion coefficients. The apparently best-match to the 
data yields a factor of 5 = 2000 faster than Bohm diffusion (with the same D oc , 6 = 1 
energy dependence). Lower exponents 6 = 0.6 produce acceptable fits with modifications of 
the model parameters within a reasonable range, whilst significantly lower values seem to 
be excluded in this scenario. We emphasize that these are illustrative models rather than 
an attempt at fully describing the system. The details depend critically on the time- and 
spatial evolution. The apparent small-scale differences between the morphology in the radio 
and GeV 7-ray bands (Figure [l] left) indicate highly non-uniform magnetic fields in Vela X. 
The true situation of escape of particles out of the ERN is likely highly complex. 



Figure [3] (left) shows the expected local electron spectrum from our model, assuming 
nominal diffusion parameters (black line, see caption for details). Several key parameters 
affect the local CR electron signature of a source: the distance, the electron injection time 
and the diffusion coefficient. We investigate the uncertainty on these parameteres by using a 
simple burst-like injection into the ISM and a spectrum adapted to match the signature of the 
full calculation (^'"^■^exp {E/5TeV) with an energy input of 6.8 x 10*^^ erg). Figure [s] (left) 
shows in red the effect of varying the time at which this burst of electrons is released into 
the ISM. It can be seen that for an injection time tinj = 3 ky, burst-like injection produces 



- 8 - 



similar results to the time-dependent injection model. Figure |3] (right) shows the effects of 
varying the distance to Vela X within the allowed errors (blue curves) (Dodson et al. 2003) 
and varying the diffusion coefficient (green curve). If source and diffusion parameters are 
changed from the nominal model it becomes difficult to keep the spectrum consistent with 
the H.E.S.S. measurements. Due to the very recent injection of a large number of high 
energy electrons, Vela X will produce a distinct feature in the local electron spectrum. Such 
a signature should be easily observable by instruments like CTA (Hermann et al. 2007), 
expected to achieve a factor > 10 increase compared to H.E.S.S. in collection area for well- 
contained/high-telescope-multiplicity events. Combined with an effective FoV solid angle 
increase of a factor > 4 and ~ 5 times longer observation time (the published H.E.S.S. data 
represents ~ 150h), the expected increase in electron statistics (with comparable systematics 
and background levels) is a factor ~ 300, extending the spectrum from ~ 4TeV up to 
~ 70TeV for an E~^ spectrum. 



4. Conclusions 

Due to its proximity and well-known distance. Vela X presents a unique opportunity to 
constrain models using excellent spectral and morphological measurements and the locally- 
measured electron/positron flux. The model presented here, of energy-dependent escape from 
Vela X, produces a clear signature in the CR electron spectrum, which could be measured 
in detail with CTA. If this scenario proves correct and Vela X proves to be typical of PWN 
of this age, the release of > TeV particles from PWN relatively early in their evolution may 
help to explain the rise in the CR positron fraction seen at high energies and at the same 
time help understand the absence of the significant population of Fermi relic-PWN that had 



been predicted before launch (see e.g. de Jager 2008). 
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Fig. 2. — Measured SEDs for the cocoon/TeV nebula (blue) and the ERN (red) compared 
to model curves. Radio data are from Alvarez et al. (2001), WMAP, ASCA and Fermi-LAT 
from Abdo et al. (2010) and TeV from Aharonian et al. (2006). The total flux of the TeV 
nebula is ~ 30% higher( Dubois et al.|2009 ), but the X-ray flux in this larger region is not well 
constrained. The magnetic fleld in the cocoon was chosen to be 4.5/iG. Injection from the 
hard X-ray nebula into this region follows a power-law with index —2 and cutoff at 70 TeV. 
The dashed line shows the ERN emission expected in the absence of escape and the solid 
curves for energy- dependent escape with 5=1000, 2000 (thick line) and 4000 times faster 
than Bohm diffusion. The current-day magnetic field is taken as 4/iG. The curves assume a 
birth-period of Pq = 22 ms and efficiency e = 0.3 (almost identical results are obtained for 
Pq = 30ms and e = 0.5) . 
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Fig. 3. — Left: Predicted electron spectrum from Vela X using the time dependent injec- 
tion model (black). Simple burst-like power-law injection spectra are shown for comparison 
assuming different injection times (all d=290 pc, and a diffusion coefficient at 10 GeV of 
5 X 10^^ cm^s^^ with energy dependence of E°'^). Right: Effects of varying distance (blue) 
or diffusion coefficient (green) for burst-like injection (at 3100 years). Data points show the 
H.E.S.S. electron measurements (Aharonian et al. 2008; Egberts et al. 2011). Shaded areas 
represent the systematic uncertainties. 



